Gap opening in graphene by shear strain 
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We exploit the concept of strain-induced band structure engineering in graphene through the calculation of 
its electronic properties under uniaxial, shear, and combined uniaxial-shear deformations. We show that by 
combining shear deformations to uniaxial strains it is possible modulate the graphene energy gap value from 
zero up to 0.9 eV. Interestingly enough, the use of a shear component allows for a gap opening at moderate 
absolute deformation, safely smaller than the graphene failure strain. 
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Graphene exhibits a number of exotic electronic properties, 
such as unconventional integer quantum Hall effect, ultrahigh 
electron mobility, electron-hole symmetry and ballistic trans- 
port even at room temperature. " Full account of these fea- 
tures is provided by the relativistic Dirac theory^ suitably de- 
veloped within the standard condensed matter formalism. A 
key feature of graphene is that its electronic density of states 
vanishes at the so-called Dirac points, where the valence and 
the conduction bands cross with a linear energy-momentum 
dispersion. Due to the hexagonal symmetry of graphene, the 
Dirac points are located at two high-symmetry points of its 
Brillouin zone. 

While many other properties of graphene are very promis- 
ing for nanoelectronics, its zero-gap semiconductor nature is 
detrimental, since it prevents the pinch off of charge current 
as requested in conventional electronic devices. Different at- 
tempts have been therefore tried in order to induce a gap, for 
instance by quantum confinement of electrons and holes in 
graphene nanoribbons^ or quantum dots. These patterning 
techniques are unfortunately affected by the edge roughness 
problem, namely: the edges are extensively damaged and 
the resulting lattice disorder can even suppress the efficient 
charge transport. The sensitivity to the edge structure has been 
demonstrated through explicit calculations of the electronic 
states in ribbons. More recently, it has been shown experi- 
mentally that a band gap as large as 0.45 eV can be opened if 
a graphene sheet is placed on an Ir(l 1 1) substrate and exposed 
to patterned hydrogen adsorption. 

Alternatively, an electronic band gap can be obtained 
by growing graphene sheets on an appropriately cho- 
sen substrate, inducing a strain field controllable by 
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Recently, it has been experimentally 



shown that by using flexible substrates a reversible and con- 
trolled strain up to ~ 18%'" can be generated with measur- 
able variations in the optical, phonon and electronic proper- 
ties of graphene. ~ This interesting result suggests that gap 
opening could be engineered by strain, rather than by pattern- 
ing. The idea has been theoretically validated by Peirera and 
Castro Neto showing that a gap is indeed generated by ap- 
plying an uniaxial strain as large as ^23%, approaching the 
graphene failure strain £y = 25%. ' This large value stands 
for the high robustness of the gapless feature of graphene un- 
der deformation. The same Authors propose an alternative 
origami technique ' ' aimed at generating local strain profiles 
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Figure 1 : (Color online) Top view of the hexagonal graphene lattice 
with its lattice vectors ai.2 = «() (^Tji'x)' where oq is the equilib- 
rium C-C distance. Axis x and y corresponds to the armchair and 
zig-zag direction, respectively. Shaded area represents the unit cell. 



by means of appropriate geometrical patterns in the substrate, 
rather than by applying strain directly to the graphene sheet. 

In this work, we further develop the above concept of strain 
engineering by showing how a combination of shear and uni- 
axial strain can be used to open a gap in a range of reversible 
and more easily accessible deformations, ranging in between 
12% and 17%. We also discuss the merging of Dirac points,'** 
which is involved into the gap opening process. 

The electronic structure of graphene has been computed for 
each deformed configuration by means of a semi-empirical 
sp^ tight-binding (TB) model, making use of the two-center 
parameterization by Xu et al. . ' Despite its semi-empirical 
character, the present TB model correctly provides the oc- 
currence of Dirac points in the band structure of graphene 
in its equilibrium geometry. Furthermore, the Xu et al. 
parametrization provides accurate scaling functions for the 
variation of the TB hopping integrals upon lattice distortions. 
This feature is instrumental for investigating gap opening in 
graphene by strain. The reliability of the present TB model 
in describing the strain-related features of graphene has been 
recently established.'^^'-" 

Graphene is an hexagonal lattice with two carbon atoms per 
unit cell and a lattice basis defined by the vectors (a 1,02), as 
shown in Fig.l, with a nearest-neighbor carbon-carbon dis- 
tance as small as aq = 1 -42 A. The in-plane elastic behavior 
of the honeycomb lattice is isotropic in the linear regime, but 
two inequivalent crystallographic directions can be neverthe- 
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Figure 2: (Color online) Brillouin zone of graphene under strain. The 
shaded areas are the corresponding irreducible part: (a) undeformed 
BZ with 6/mmm hexagonal symmetries; (b) BZ deformed by uniax- 
ial strain with mmm rhombic symmetry; (c) BZ deformed by shear 
strain with 2/m monoclinic symmetry. 



less defined: the so-called armchair and zig-zag directions, 
shown in Fig.l as x and y axis, respectively. According to 
the Cauchy-Born rule, when straining a graphene sample its 
lattice vectors are affected accordingly, as well as the asso- 
ciated reciprocal vectors ^^1,^2^ . The deformed vectors are 

given by a'- = {EijUj + dijaj), where e = {£;;} is the strain 
tensor describing the deformation and i,j = x,y. The condi- 
tion a'l^-b'i— 2n5iii (where ^, / = 1 , 2) allows one to obtain the 
deformed reciprocal lattice vectors. 

The following in-plane deformations have been applied to 
the equilibrium honeycomb lattice under plane-strain border 
conditions: (i) an uniaxial deformation ^ along the armchair 

direction, corresponding to a strain tensor Sj"'^^ = ^Sixdj^, (ii) 
an uniaxial deformation ^ along the zig-zag direction, cor- 
responding to a strain tensor E/j^'^ = i^5,y5,y; (iii) an hydro- 
static planar deformation ^, corresponding to the strain tensor 
' — C^^ ; (i^) ^ shear deformation ^, corresponding to an 

in-plain strain tensor e^-^^' = ^ {Six5jy + 5iy5jx)- 

In order to extend the reliability of the present model to 
electronic features under strain, our results about the effects 
of hydrostatic and uniaxial deformations on the band structure 
are at first compared with previous data available in literature. 
For graphene under in-plane hydrostatic deformation with 
^ < 15%, both in compression and in traction, we have cal- 
culated the band electronic structure and the density of states. 
Since the hydrostatic strain does not change the Dg/, (6 / mmm) 
symmetry of the hexagonal lattice (Fig. 2a), we only observe 
the variation of the pseudogaps at F and M points, while the 
location of the Dirac points is clamped at the K point. In par- 
ticular, the pseudogap at M decreases almost linearly from 6 
eV (for C = -15%) to 1.8 eV (for C = +15%). We remind 
that its value for the unstrained configuration is 2.2 eV. These 
results are in quantitative good agreement with Ref.2 1 . Any 
other non-hydrostatic deformation lowers the symmetry of the 
graphene lattice. When an uniaxial strain is applied, all the 6- 
and 3-fold rotational symmetries are lost: a transition from the 
hexagonal Dg/, (6 / mmm) to the rhombic D2/, (mmm) symmetry 
is observed (Fig. 2b). The irreducible part of the first Bril- 
louin zone (BZ) is also affected by such deformations, since 
its original triangular shape (Fig. 2a) is varied to the polygo- 
nal form represented in Fig. 2b. The top of the valence band 
and the bottom of the conduction band are shown in Fig. 3 for 



(a) undeformed graphene 
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(b) uniaxial (armchair) strain with ^ = 15% 




(c) uniaxial (zig-zag) strain with f = 15% 
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Figure 3: (Color online) Top of the valence band (red, marked as 
VB) and bottom of the conduction band (green, marked as CB) of 
graphene under uniaxial strain. Panel a: band structure of the unde- 
formed lattice. Panel b and c: band structure under uniaxial strain 
along the armchair and the zig-zag direction, respectively. Symbols 
connect the high-symmetry points of the BZ (bottom shaded area) to 
the energy of the corresponding electronic states. 



the undeformed configuration (panel a), as well as under uni- 
axial deformation (panels b and c, corresponding to a strain 
^ — 15% along the armchair direction and in the zig-zag di- 
rection respectively). The main effect of strain is the opening 
of a pseudo-gap at K and K' . Accordingly, the Dirac points are 
no more located at such high-symmetry points; rather, they 
drift away within the BZ, either for deformations along arm- 
chair direction or along zig-zag one. Once again this impor- 
tant qualitative feature is in good agreement with Ref. 14. 

Let us now consider the case of an in-plane shear deforma- 
tion, described by the following shear strain 



where ^ is the strain parameter. Such a deformation modifies 
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Figure 4: (Color online) Top of the valence band (red, marked as 
VB) and bottom of the conduction band (green, marked as CB) of 
graphene under pure shear strain with ^ = 20%. Symbols connect the 
high-symmetry points of the BZ (bottom shaded area) to the energy 
of the corresponding electronic states. 



the original reciprocal lattice vectors h\ and into 
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By applying the shear strain given in Eq. 1 to the graphene 
lattice, its symmetry class is further lowered to monoclinic. 
The corresponding symmetry group is 2 jm. Because of this 
change in symmetry, the irreducible part of the BZ is affected 
accordingly as shown in Fig. 2c, which has K', K, K", K' as 
corners. In the undeformed lattice, a Dirac point is located 
at each of these corners. The scenario under shear strain is 
quite different from the case of uniaxial deformations: at the 
comparatively small strain ~ 16%, a gap is indeed opened. 
The rise of a gap in the electronic band structure under shear 
is due to a peculiar process that involves the merging-- of two 
Dirac points, namely D' and D", which move away from the 
corners K' and /T" and approach each other inside the BZ. 

Is important to remark that the merging of the two inequiv- 
alent Dirac points and the opening of a gap, appear for a shear 
strain value ^ ~ 16% which is lower than in the case of zig-zag 
uniaxial deformation. The gap increases up to a maximum 
value of 0.72 eV for shear strain parameter of ~ 20%, as 
shown in Fig. 4. We conclude that shear strain seems a likely 
candidate to achieve gap opening in graphene for a deforma- 
tion far enough from failure strain and, therefore, achievable 
with no danger for the overall mechanical stability of the two- 
dimensional sheet. 

Gap opening is predicted by the present TB calculation to 
occur at an even smaller strain parameter i^, provided that a 
combination of shear and uniaxial strain is considered. By 
adding an uniaxial component to shear we generate a strain 
tensor of the form 
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for which the symmetry class of the lattice is not changed with 
respect to the pure shear case. 



Nevertheless, uniaxial deformations along the armchair or 
zig-zag direction are found to dissimilarly affect the band 
structure of graphene. Only in the last case we have observed 
the merging of the Dirac points already at ^ ~ 12%. The main 
features of the transition is the same as described before. The 
energy gap grows up to a maximum value of 0.95 eV (when 
the strain parameter achieves a value of ^ ~ 17%), reducing 
again to zero at ^ ~ 20% due to the steady decrease of the 
direct gap at F . 

In order to quantitatively describe the evolution of the gap 
opening as function of the applied strain, the density of states 
(DOS) has been calculated by a two-dimentional 75 x 150 x 1 
regular ^-point mesh of the (deformed) BZ. As shown in Fig. 5, 
for a strain value less then 15% (panel a) or 11% (panel b), 
the DOS depends linearly on energy close to the Fermi level, 
showing a slope increasing with the strain. The two charac- 
teristic Van Hove singularities into the DOS move closer the 
Fermi energy and disappear into abrupt gap-edges as soon as 
the gap is open. After the annihilation of the Dirac points, the 
DOS shows a ~ \fE behavior 

We conclude by remarking that the two strain contributions 
(i.e. uniaxial and shear) could be combined in different ways 
so as to modulate the energy gap value. In Fig. 6, the elec- 
tronic band structures of graphene under different combina- 
tions of shear and uniaxial strain are compared, keeping the 
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Figure 5: (Color online) Density of states around the Fermi level (set 
conventionally at eV) as function of the strain parameter C, . Panel 
(a); graphene under pure shear deformation. Panel (b): graphene 
under combined shear and uniaxial deformation (along the armchair 
direction). The maximum value of the energy gap is observed for a 
strain parameter as large as ~ 20% and [, ~ 17% respectively. 
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(a) (shear + zig-zag uniaxial) deformation 



(b) (shear + armchair uniaxial) deformation 
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Figure 6: (Color online) Top of the valence band (red, marked as VB) and bottom of the conduction band (green, marked as CB) of graphene 
under combined shear and uniaxial strain with ^ = 15%. The uniaxial component of the strain is applied along the zig-zag (panel a) and 
armchair (panel b) direction. Symbols coimect the high-symmetry points of the BZ (bottom shaded area) to the energy of the corresponding 
electronic states. 



same value of the strain parameter ^ = 15%. While the com- 
bination of shear with uniaxial armchair shows a sizable en- 
ergy gap of about 0.6 eV, the combination of shear with uni- 
axial zig-zag is associated to a gapless band structure. 

In conclusion, we have shown how the opening gap in the 
electronic spectrum of graphene could be achieved by apply- 
ing deformations with a nonzero shear component, rather than 
a simple uniaxial deformation. Energy gaps have been found 
to vary in between 0.6 eV and 0.9 eV for a strain value far 



enough from failure. In particular, we have shown that the 
most effective way to control the gap opening is to combine a 
shear and an armchair uniaxial deformation. 

We acknowledge computational support by COSMOLAB 
(Cagliari, Italy) and CASPUR (Rome, Italy). Discussions 
with G. Fadda are gratefully acknowledged. One of us (L.C.) 
acknowledges partial financial support by the MATHMAT 
project (Universita di Padova, Italy). 



* E-mail me at:luciano. colombo@dsf.unica.it 
' K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 1 Y. Zhang, 
S. V. Dubonos, I. V. Grigorieva, A. A. Firsov, Science 306, 666 

(2004) . 

^ A. H. Castro Neto, F.Guinea, N. M. R. Peres, K. S. Novoselov, 
and A. K. Geim, Rev. Mod. Phys. 81, 109 (2009). 

^ K. S. Novoselov, Z. Jiang, Y. Zhang, S. V. Morozov, H. L. 
Stormer, U. Zeitler, J. C. Maan, G. S. Boebinger, P. Kim, A. K. 
Geim, Science 315, 1379 (2007). 

V.R Gusynin and S.G. Sharapov, Phys. Rev. Lett. 95, 146801 

(2005) . 

^ M. Y. Hani, B. Ozyilmaz, Y. Zhang, and P. Kim, Phys. Rev. Lett. 

98, 206805 (2007). 
^ F. Sols, F. Guinea and A. H. Castro Neto, Phys. Rev. Lett. 99, 

166803 (2007). 

' E. R. Mucciolo, A. H. Castro Neto and C. H. Lewenkopf, Phys. 

Rev. B 79,075407 (2009). 
^ K. Nakada, Fujita, G. Dresselhaus and M.S. Dresselhaus, Phys. 

Rev. B 54, 17954 (1996). 
^ Z. H. Ni, H. M. Wang, Y. Ma, J. Kasim, Y. H. Wu and Z. X. Shen, 

ACS Nano 2, 1033 (2008). 

R. Balog, B. J0rgensen, L. Nilsson, M. Andersen, E. Rienks, M. 
Bianchi, M. Fanetti, E. Lsegsgaard, A. Baraldi, S. Lizzit, Z. Slji- 
vancanin, F. Besenbacher, B. Hammer, T. G. Pedersen, P. Hof- 
mann and L. Hornekaer, Nature 9, 315 (2010). 

" P Shemella and S.K. Nayak, Appl. Phys. Lett. 94, 032101 (2009). 

12 K. S. Kim, Y. Zhao, H. Jang, S. Y. Lee, J. M. Kim, K.S. Kim, J. H. 
Ahn, P Kim, J. Y. Choi and B. H. Hong, Nature 457, 706 (2009) 
Z. H. Ni, T. Yu, Y. H. Lu, Y. Y. Wang, Y. P Feng and Z. X. Shen, 



ACS Nano 2, 2301 (2008). 
1"* V. M. Pereira and A. H. Castro Neto, Phys. Rev. B 80, 045401 
(2009). 

15 C. Lee, X. Wei, J. W. Kysar and J. Hone, Science 321, 385 (2008). 
E. Cadelano, P. L. Palla, S. Giordano and L. Colombo, Phys. Rev. 
Lett. 102, 235502 (2009). 

V. M. Pereira and A. H. Castro Neto, Phys. Rev. Lett. 103, 046801 
(2009). 

G. Montambaux, F. Piechon, J.-N. Fuchs and O. Goerbig, Phys. 

Rev. B 80, 153412 (2009). 
■■^ C. H. Xu, C. Z. Wang, C. T. Chan and K. M. Ho, J. Phys.: 

Condens. Matter 4, 6047 (1992). 
2" E. Cadelano, S. Giordano and L. Colombo, Phys. Rev. B 81, 

144105 (2010). 

2' G. Gui, J. Li and J. Zhong, Phys. Rev. 78, 075435 (2008) 
22 By merging we mean the following: at a critical strain ^ ~ 
15.95% the Dirac points are so close that they annihilate in a sin- 
gle hybrid Dirac cone, which shows a peculiar energy-momentum 
dispersion: it is linear (quadratic) along (perpendicular) to the di- 
rection joining the two Dirac points. At = 16.0%, a gap as small 
as 0.05 eV is eventually opened. More details about the motion of 
Dirac points in two-dimensional crystals under uniaxial stress are 
reported in Ref. 18,23,24. 

P. Dietl, F. Piechon and G. Montambaux, Phys. Rev. Lett. 100, 
236405 (2008). 

2"* Y. Hasegawa, R. Konno, H. Nakano and M. Kohmoto , Phys. Rev. 
B 74, 033413 (2006). 



